A relativistic jetted outburst from a massive black hole 
fed by a tidally disrupted star 



Joshua S. Bloom 1 *, Dimitrios Giannios 2 , Brian D. Metzger 2,3 , 
S. Bradley Cenko 1 , Daniel A. Perley 1 , Nathaniel R. Butler 13 , 

Nial R. Tanvir 4 , Andrew J. Levan 5 , Paul T. O' Brien 4 , 
Linda E. Strubbe 1 ' 6 , Fabio De Colle 7 , Enrico Ramirez-Ruiz 7 , 

William H. Lee 8 , Sergei Nayakshin 4 , Eliot Quataert 1,6 , 
Andrew R. King 4 , Antonino Cucchiara 1 ' 9 , James Guillochon 7 , 
Geoffrey C. Bower 101 , Andrew S. Fruchter 11 , Adam N. Morgan 1 , 
Alexander J. van der Horst 12 , 

department of Astronomy, University of California, Berkeley, CA 94720-3411, USA 
2 Department of Astrophysical Sciences, Peyton Hall, Princeton University, Princeton, NJ 08544, USA 

3 NASA Einstein Fellow 

department of Physics and Astronomy, University of Leicester, University Road, Leicester LEI 7RH, UK 
department of Physics, University of Warwick, Coventry CV4 7AL, UK 
6 Theoretical Astrophysics Center, University of California, Berkeley, CA 94720, USA 
7 Astronomy and Astrophysics Department, University of California, Santa Cruz, CA 95064, USA 

8 Instituto de Astronomi'a, Universidad National Autonoma de Mexico, Apdo. Postal 70-264, Cd. Universitaria, Mexico DF 04510 
Computational Cosmology Center, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA 
10 Radio Astronomy Laboratory, University of California, Berkeley, 601 Campbell Hall 341 1, Berkeley, CA 94720, USA 

n Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 
^Universities Space Research Association, NSSTC, 320 Sparkman Drive, Huntsville, AL 35805, USA 
*To whom correspondence should be addressed; E-mail: jbloom@astro.berkeley.edu. 

While gas accretion onto some massive black holes (MBHs) at the centers of 
galaxies actively powers luminous emission, the vast majority of MBHs are 
considered dormant. Occasionally, a star passing too near a MBH is torn apart 
by gravitational forces, leading to a bright panchromatic tidal disruption flare 
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(TDF) (1-4). While the high-energy transient Swift J164449.3+573451 ("Sw 
1644+57") initially displayed none of the theoretically anticipated (nor previ- 
ously observed) TDF characteristics, we show that the observations (5) sug- 
gest a sudden accretion event onto a central MBH of mass ~ 10 6 — 10 7 solar 
masses. We find evidence for a mildly relativistic outflow, jet collimation, and 
a spectrum characterized by synchrotron and inverse Compton processes; this 
leads to a natural analogy of Sw 1644+57 with a smaller-scale blazar. The phe- 
nomenologically novel Sw 1644+57 thus connects the study of TDFs and active 
galaxies, opening a new vista on disk-jet interactions in BHs and magnetic field 
generation and transport in accretion systems. 

While variability is common to all active galactic nuclei (AGN) — fundamentally tied to the 
unsteady accretion flow of gas towards the central MBH — the timescale for active MBHs to dra- 
matically change accretion rates (leading the source to, for example, turn "off"), is much longer 
than a human lifetime. The most variable AGN are a subclass called blazars, with typical masses 
M B h ~ 10 8 — 1O 9 M (M Q is the mass of the Sun), originally found to be radio and optically 
bright but with luminosities dominated by X rays and gamma rays. Significant changes in the 
apparent luminosity over minutes- to hour-long timescales are thought to be predominately due 
to Doppler-beamed emitting regions within a jetted outflow moving relativistically toward the 
observer (6, 7). The typical Lorentz factor inferred is Tj = (1 — /3 2 )" 1 / 2 10, with velocity 
v = (5c of the jetted outflow. The high-energy emission is thought to be due to inverse Compton 
upscattering of the accretion disk photons, photons from within the jet itself, and/or photons 
from structures external to the accretion disk (8, 9). 

Inactive MBHs can suddenly "turn on" while being fed by temporary mass accretion es- 
tablished following the tidal disruption of a passing star (1). If a star of mass M* and radius 

passes within the disruption radius r d ps /^(Mbh/M*) 1 / 3 ps 5Mf 2 ^r s (with M B h = 

2 



1O 7 M 7 M and r s = 2GM BR /c 2 the Schwarzschild radius of the BH, M* = M , = R e ), 
then a mass of ~ M*/2 will accrete onto the MBH with a peak accretion rate on a timescale 
of weeks (2). Importantly, the accretion rate for typical scenarios with a M 7 BH can be super- 
Eddington for months (4, 10). Candidate TDFs have been observed at X-ray, ultraviolet, and 
optical wavebands (11-14), but the rates are very low [~ 1CT 5 yr gal -1 ; (3)] and the connection 
between the observed light curves and spectra to theoretical expectations has been tenuous. 

Recently, Giannios & Metzger (15) suggested that a ~ 1 yr radio event could follow a 
TDF arising from a jetted relativistic outflow as it interacted with (and was slowed down by) 
the external ambient medium, akin to the afterglow from external shocks following gamma-ray 
bursts (16). The supposition was that the observer viewed the event off-axis from the relativistic 
jet. Just what would be seen if instead the jet were pointed nearly towards the observer — as in 
the geometry inferred for blazars — was not considered. 

Sw 1644+57: a relativistic outflow generated by an accreting 
massive black hole 

Sw 1644+57 was initially discovered as a long-duration gamma-ray burst (GRB 110328A) by 
the Swift satellite (17) at a time t = 201 1 March 28 12:57:45 UT but, given the longevity and 
flaring of the X-ray afterglow, it was quickly realized that the high-energy emission was unlike 
that associated with any GRB previously; in a companion paper, Levan et al. (5) (who find 
redshift of the host galaxy of Sw 1644+57 to be z =0.3543) describe in detail the differences 
between this event and the GRB population observed by Swift. Based on the data available 
in the first two days following the event, we (18) first suggested^ a possible analogy of Sw 
1644+57 with a scaled-down version of a blazar impulsively fed by ~ 1M of stellar mass. 

There are several lines of evidence to suggest an accreting MBH origin. First, the astro- 

1 A similarly detailed sketch was also posted (19) later. 
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metric coincidence of the X-ray, optical, infrared, and radio transient with the light-centroid of 
the putative host galaxy is strongly indicative of a positional connection to an MBf^] Second, 
the compact nature of the associated radio source suggests an emission region no larger than 5 
parsec which disfavors other possible extended emission sources (see SOM). Last, the observed 
correlation between the X-ray flux and spectral hardness (see SOM) is similar to that observed 
in blazars (21, 22). 

Accepting the accreting MBH hypothesis, we now examine constraints on the BH mass 
and the accretion characteristics. We infer a minimum host-frame variability timescale of 
£var,min ~ 78 sec from the X-ray light curve (SOM). By requiring that t var ,min exceed the light- 
crossing time of r s , we derive an upper limit on the MBH mass Mbh ^8x 1O 6 M . Note that, 
like in gamma-ray burst light curves, even in the presence of relativistic motion (see below), 
the observed variability should track that of the energy injection timescales from the central 
engine. Irrespective of the timing argument, we can place approximate upper limits on the mass 
of the central BH if we assume the entirety of the galaxy mass [fewxlO 9 M ; (5)] and light 
[fewxlO 9 L ; (5)] arise from the host bulge (i.e., not in the disk) and apply the bulge mass- 
BH mass and the bulge luminosity-BH mass correlations (23-25). All such analyses suggest 
.Mbh ^ 1O 7 M , securely under the limit (fewxlO 8 M ) required for the tidal disruption of a 
solar-mass star to occur outside the event horizon of the MBH. 

If the emission is isotropic, the average X-ray luminosity of the outburst (SOM), L x ~ 10 47 
erg s _1 , corresponds to the Eddington luminosity of a ~ 1O 9 M BH, incompatible with the 
upper limit derived from variability. If the source is relativistically beamed (see SOM), with 
beaming factor f b — (1 — cos6>j) < 1, the beaming-corrected luminosity fbL x ~ 10 45 erg s _1 
becomes consistent with the Eddington luminosity of a ~ 1O 7 M SMBH if 9j = 1 /Tj m 0.1, as 

2 Indeed we predicted the precise astrometric coincidence in (18) before the connection was solidified (20). 
However, within the uncertainties from Hubble imaging [^300 pc; (5)], the central stellar and gas density could 
be high enough to allow other progenitors, such as supernovae. 
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inferred in blazars (we show below that this value of Tj is also consistent with the inferred rate 
of Sw 1644+57-like events). We can also infer the presence of relativistic outflow (see SOM) 
by requiring that the true brightness temperature of the radio transient be less than the inverse 
Compton catastrophe temperature 10 12 K. Those constraints require a mean Tj ^ 1.9 from t to 
the time of our VLBI observations (5). Separately, we can use the observed variability of the 
radio counterpart to place constraints on the source size, finding Tj ^ 10. 

A tidal disruption origin for the accretion mass 

If the source had been active in the distant past, we would expect to observe extended ra- 
dio emission (e.g., jets or other emission knots) in our VLBI imaging. Since this is not seen 
and archival searches spanning two decades have yielded no evidence for prior AGN activity 
from radio to gamma-ray wavebands (see SOM), the evidence thus suggests that a M B h = 
10 6 — 1O 7 M BH underwent a dramatic turn on to near-Eddington accretion rates, launching 
an energetic, relativistic outflow in the process. This rapid increase in the accretion rate cannot 
result from gas entering the sphere of influence (soi) of the MBH, since this would require a 
timescale ~ R SO i/<y ^ 10 4 yr to appreciably alter the accretion rate near the horizon, where 
-R soi ~ 1 pc is the radius of the sphere of influence and a ~ 100 km s _1 is a typical bulge veloc- 
ity dispersion. We suggest instead that a TDF provides a natural explanation for Sw 1644+57. 

We again turn to the energetics and timescales to provide guidance on the TDF hypothesis. 
The observed X-ray fluence Sx suggests an energy release of E x = 9.2 x lO 52 /^ erg for the 
first 10 6 s. Assuming the source continues to radiate at a sustained level for ~1 month, that the 
energy released in the XRT band is about 1/3 of the bolometric energy (following from Fig. 1) 
and adopting f b = 5x 10~ 3 , the total energy release from the jet amounts to 0.3% of the mass- 
energy accreted if M* = M Q . Given a typical accretion efficiency of e B H = E av /m acc c 2 = 0.1, 
the jet need radiate only about l/30th of the available energy E av . The duration of the X-ray 
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light curve and the requisite accretion rate are also compatible with the several-day fallback 
timescale (see SOM). 

The broadband SED of Sw 1644+57 shown in Fig. 1 displays two peaks, at far infrared and 
X-ray /gamma-ray wavebands, respectively. The overall spectral shape is reminiscent of blazars, 
for which the peaks at low and high energies are typically modeled as synchrotron and Inverse 
Compton (IC) emission, respectively. The X-ray emission shows both a bright/flaring and a 
dim/slower-varying ("quiescent") state. Under the TDF hypothesis, what could account for the 
observed spectrum and temporal behavior? 

• Single Component Synchrotron with Dust Extinction: In low-luminosity BL Lac ob- 
jects, the vF v synchrotron spectral peak may occur at energies as high as hard X-rays. 
Thus, one possibility is that the entire emission from radio to X-rays is part of a single 
non-thermal synchrotron spectrum originating from shocked relativistic electrons. In this 
scenario, the suppressed optical emission and red IR colors of the transient could result 
from dust extinction with Ay > 10 mag. We conclude that although a single extinguished 
synchrotron spectrum cannot be ruled out, the large required extinction may disfavor this 
interpretation (see Fig. S3). Furthermore, although a synchrotron origin is still likely for 
at least the radio emission, there is evidence that the radio and X-ray emitting regions 
may not be coincident (see below). 

• Two- Component Blazar Emission: The FIR and hard X-ray peaks may, instead, repre- 
sent distinct spectral components, corresponding to synchrotron and IC emission, respec- 
tively, as in blazars. Figure 1 shows that the vF v luminosity of the low energy peak is 
~ 1 — 2 orders of magnitude weaker than the high energy peak. This extreme ratio, and 
the relatively low frequency of the synchrotron peak, are both compatible with Eddington- 
accreting blazar emission (9). We consider various origins for the seed photons for IC in 
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the SOM and show two example model SEDs in Figure 1. 

• Forward Shock Emission from Jet-ISM Interaction: Although the above models gen- 
erally assume that the low- and high-energy spectral components are directly related, evi- 
dence suggests that they may originate from distinct radii, at least during the X-ray flaring 
state. While the rapid variability of the X-ray emission strongly indicates an "internal" 
origin (26), the radio-IR emission varies more smoothly and could instead result at larger 
radii from the interaction of the jet with the surrounding interstellar medium (see SOM). 
If no AGN activity occurred prior to the recent onset of emission, the jet must burrow 
its way through the gas in the nuclear region^] Due to its fast motion, the newly-formed 
jet drives a shock into the external gas (forward shock), while simultaneously a reverse 
shock slows it down. Particles accelerated at these shocks may power synchrotron af- 
terglow emission beginning simultaneously when the jet forms, yet lasting long after the 
internal emission has faded. This model, the geometry of which is depicted in Figure 2, 
appears to best accommodate the data, and makes specific predictions for the long-term 
fading of the radio and IR transient (see SOM). 

Conclusions and Predictions 

No rising UV-optical transient nor slowly evolving thermal X-ray component has been seen to 
date from Sw 1644+57, in contrast with the nominal expectations of TDFs. However, if Sw 
1644+57 was obscured by dust, then UV-optical suppression of the transient would be natural. 
And if we understand the thermal X-ray emission as being outshone by the jetted emission in the 
first week, then the TDF hypothesis would naturally lead to the prediction of the emergence of 
the thermal component on timescales of weeks to months. Whether and when thermal emission 

3 This situation is not encountered in normal (long-lived) blazars because a large ;> kpc scale cavity has been 
carved by the preceding outflow. 
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will be observable hinges on the degree of dust extinction, its brightness relative to the host 
bulge, and how rapidly the jet emission fades. 

If the TDF hypothesis is correct, Sw 1644+57 will fade over the coming year and will 
not repeat. If our interpretation about the relativistic flow and spectral origin is correct, then 
we would expect the transient emission to be polarized at a (low) level similar to that seen in 
gamma-ray burst afterglows (as opposed to blazars^]). Moreover we expect to see evidence for 
superluminal motion of the radio source as seen in VLBI monitoring over the next few months; 
the source itself may become resolved on timescales of a few months if it remains bright enough 
to detect at radio wavebands. 

Adopting a beaming fraction fa ^ 10~ 2 consistent with that inferred from Sw 1644+57 
(SOM), we conclude that for every on-axis event, there will be 1/ f b ^ 10 2 events pointed away 
from our line of sight. Since Swift has triggered on one event in ~ 6 years of monitoring, the 
total inferred limit on the rate of TDFs accompanied by relativistic ejecta is ^ 10 yr 1 out to 
a similar distance. Although the majority of such events will not produce prompt high-energy 
emission, bright radio emission is predicted once the ejecta decelerates to non-relativistic speeds 
on a timescale ~ 1 year (15). The predicted peak flux is sufficiently high (~ 0.1 — 1 mJy at 
several GHz frequencies and redshifts similar to Sw 1644+57) that ~ 10 — 100 may be detected 
per year by upcoming radio transient surveys. Continued long-term monitoring of the radio 
afterglow from Sw 1644+57 will allow for calibration of the off-axis models. 

We conclude with two broader ramifications that stem from our understanding of the origin 
of Sw 1644+57. The emerging jet from the tidal disruption event appears to be powerful enough 
to accelerate cosmic rays up to ~ 10 20 eV, i.e., the highest observed energies (28). The observed 

rate of jets associated with the tidal disruption of a star, R ~ 10 _11 Mpc~ 3 yr _1 , and the energy 

4 Here the departure from the blazar analogy is worth noting in that the physics of the radio emission is likely to 
be different in this case: we have suggested that the emission is originating from the shocked surrounding material 
(forward shock), not the shocked jet as in normal blazars, which could contain large scale fields. Even so, only 
10% of flat-spectrum radio quasars and BL Lac objects have polarization larger than our VLBI limits (27). 
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released per event of E x ~ 10 53 erg, however, imply an energy injection rate of -Stdf ~ 
10 42 erg Mpc _3 yr _1 . Despite the large uncertainty, this rate is significantly smaller than the 
injection rate E hl j ~ 10 44 erg Mpc _3 yr _1 required to explain the observed flux of cosmic rays 
of energy > 10 19 eV (29). This conclusion is, however, subject to uncertainties associated with 
the radiative efficiency of the jet. 

There is much evidence that AGN jets are accelerated by magnetohydrodynamic, rather 
than hydrodynamic, forces (30). A key unsolved question is whether the large-scale magnetic 
field necessary to power the jet is advected in with the flow (31), or whether it is generated 
locally in the disk by instabilities or dynamo action. If the jet is launched from a radius R m , the 
magnetic field strength at its base (B) is related to the jet luminosity by L } ~ nRf n c x (B 2 /4ir). 
If we assume L } ~ 10 45 ergs s _1 , similar to the Eddington limit for a ~ 1O 7 M MBH (as 
appears necessary to explain the bright non-thermal emission), the required field strength is 
B ~ 10 5 G for R m ~ 1.5r s . This field is much higher than the average field strengths of typical 
main sequence stars (< 10 3 G). The stellar field is further diluted due to flux freezing by a factor 
~ (R*/R- ni ) 2 as matter falls into the BH, where R± ~ R Q is the stellar radius prior to disruption. 
Hence we conclude that the large-scale field responsible for launching the jet associated with 
Sw 1644+57 must have been generated in situ. Placing similar constraints has not previously 
been possible in the context of normal AGN or X-ray binary disks, due to the much larger ratio 
between the outer and inner disk radii in these systems. 
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Observed frequency V (Hz) 



Fig. 1: Multiwavelength spectral energy distribution of Sw 1644+57 at to + 2.9 day. Our 
radio-through-UV measurements are represented by filled circles, with data from the pub- 
lished circulars (5) represented by open circles; the uncertain relative contributions of the 
host galaxy and the optical/infrared (IR) transient result in very large uncertainties for the 
J and z photometric data points. X-ray and soft gamma-ray points from the Swift XRT 
and BAT (uncorrected for host-galaxy absorption) are shown as black crosses, and the 
Fermi/LAT gamma- ray upper-limit (32) is shown at the far right. The 90% uncertainty 
region of a power-law fit to the XRT data (with N H absorption removed) is represented by 
the blue bow-tie. The inset at lower right shows the same data zoomed in on the optical- 
NIR window. These observations are overplotted with two different multi-component 
models for the SED, exploring different emission mechanisms and radii (see also Fig. 2). 
The orange curve shows a model with synchrotron, synchrotron self-Compton, and ex- 
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ternal Compton contributions. In this model, the radio and IR emission are produced by 
synchrotron radiation from an extended source, while the X-ray emission is dominated 
by the Compton scattering of external photons from the accretion disk (for illustrative 
purposes, we assume a 1O 6 M MBH). The purple curve shows a model in which the IR 
emission originates from a compact source of synchrotron emission, as may be required 
by the IR variability (~ 4 x 10 14 cm). As in the orange model, the X-ray emission is 
dominated by external Compton scattering, while the peak at high energies results from 
synchrotron self- Compton emission. An additional synchrotron component from a mildly 
relativistic blast- wave afterglow at larger radius is invoked to explain the bright radio and 
millimeter fluxes. Both models require moderate extinction (A v ~ 3 — 5 mag). Additional 
synchrotron models are shown in Figure S3. The model SEDs in this Figure and in Figure 
S3 were generated using the computer code from Ref. (33). 
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Sw 1 644+57 




Fig 2. Schematic representation of the geometry and emission regions for Sw 1644+57. A 
star is disrupted at distance r d from a black hole of mass M B h with Schwarzschild radius 
r s . Half of the mass of the star escapes on unbound orbits while the other half remains 
bound. Shocked, circularized fallback mass sets up a temporary accretion disk with in- 
ner radius 3r s (for a non-spinning BH). A two-sided jet is powered starting at the time 
of accretion and plows through the interstellar region surrounding the BH at a Lorentz 
factor Tj. At some later time, the jet has reached a distance R s where the forward shock 
radiates the observed radio and infrared light. Emission from the accretion disk is Comp- 
ton upscattered giving rise to the observed X-rays. Different viewing angles (whether the 
observer is inside 0j 1/Tj or not) determines what sort of phenomena is observed. An 
analogy with blazars and AGN for more massive BHs is given. 
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1 Supplemental Online Material 



"A relativistic jetted outburst from a massive black hole fed by a tidally dis- 
rupted star," by J. S. Bloom et al. 

In what follows, we assume cosmological parameters H = 73 km s -1 Mpc -1 , 
tt A = 0.73 and VL m = 0.27. At a redshift of z = 0.3534 (1), the luminosity 
distance is 1814.8 Mpc and the angular diameter distance cLa = 990.8 Mpc (with 
1 arcsec corresponding to 4.803 kpc in projection). 

1.1 Timing and Energetics Analysis of the X-ray Light Curve 

Figure SI shows the X-ray light curve in the first 16.6 days from the Swift/BAT 
trigger at time to- The total on-source integration time of the XRT was ~ 3 d. 
From the XRT light curve we can infer a minimum timescale of variability over 
which the flux changes by a factor of 2. We find this minimum observed t var ,mm = 
105 sec occurs as a rise between to + 1 1 1499 and t + 1 1 1604 sec. The most rapid 
observed decline (t var = 188 s) occurred much later, at a time t + 1418381 and 
t + 1418569 sec, as the XRT light curve began at major secular downturn. From 
a periodogram analysis, we also (2, 3) find no significant frequencies in the XRT 
light curve that are not associated with the orbital period of the spacecraft. 

As in accreting BH systems the X-ray hardness tracks the X-ray flux of Sw 
1644+57 (4, 5), such that the best-fit powerlaw index increases with decreasing 
flux (1). From structure function analysis (Fig. S2), we measure a red-noise power 
spectrum, characterized by a powerlaw slope with frequency a = —1.6 ± 0.2. This 
is consistent with the typical value [a = —2.0 ± 0.5; (<5)] at X-ray wavebands for 
blazars. Such a red noise spectrum is also typical of low-frequency (/ < 1 Hz) 
X-ray variability in BH and neutron star binaries (7). We thus interpret the power 
density spectrum of the X-rays as indicative of accretion around a central compact 
mass. 

The peak luminosity observed in the X-rays is Lx, pca k ~ 3.0 x 10 48 erg s' 1 and 
the average over the first 1.2 x 10 6 seconds is Lx ~ 9.2 x 10 46 erg s _1 . The fluence 
within 0.3 keV to 10 keV — which is ^20 times that from longer wavelengths (see 
Fig. 1 of the main text and Fig. S3)— from t to t + 1.2 x 10 6 s is S x = 3 x 10~ 4 
erg cm -2 . 
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1.2 Archival Search for AGN activity 

One possible explanation for the observed activity is that Sw 1644+57 represents 
an unusual AGN that went undetected up to this point. However, archival Swift 
observations place constraints on existing bright gamma-ray and X-ray emission 
dating back to Feb. 2005 (8), with typical one-day average luminosities L x ,p rc <^ 
1.7 x 10 44 erg s _1 , a factor of about 10 3 less than the average luminosity of Sw 
1644+57. The source was also not detected in the ROSAT All-Sky Survey. The 
upper limit of 1.7 x 10 44 erg s _1 implies no similar activity back to 1990 (Pilar 
Esquej, private communication). Furthermore, the FIRST survey places a 1.4 
GHz (3 a) upper limit of <; 0.3 mJy on quiescent radio emission (9). We know 
of no pre-event spectrum of the host galaxy, but post-event spectroscopy analysis 
showed line-ratio characteristics of a star-forming galaxy without any indication 
of AGN activity (1). 

1.3 Inference of Relativistic Outflow & Beaming 

High-Energy Constraints: Beaming of the observed X-ray emission (10) could 
be due either to purely geometric effects associated with the accretion disk and its 
outflows (11, 12) or due to relativistic effects in a collimated jet. Although both 
possibilities are plausible a priori, the presence of bright non-thermal radio emis- 
sion strongly suggests the presence of some relativistic material. Geometrically- 
beamed emission could be observed, even if the observer is not directly along the 
line of sight of the relativistic jet. However, the fact that the beaming fraction 
f b ~ 9f /2 ~ 5 x 1CT 3 of typical blazar jets (with 0j ~ 0.1, the half opening angle; 
Fig. 2 main text) is similar to the beaming estimated from the predicted TDE rate 
(§1.5) suggest that this possibility is not essential. The beaming-corrected lumi- 
nosity in this case fbLx ~ 10 45 erg s _1 is also then consistent with the Eddington 
luminosity of a ~ 1O 7 M MBH. 

Another possible explanation for very bright X-ray emission is that the accre- 
tion rate, and the resultant jet power, is above the Eddington rate. Although super- 
Eddington accretion is indeed predicted at early times in TDE associated with 
low SMBH masses and super-Eddington jets may have empirical support (13), 
completely unbeamed emission is inconsistent with the rate estimates (§1.5). 
Radio Constraints: We can constrain the outflow velocity of Sw 1644+57 by 
making use of the well-known brightness temperature limits on synchrotron radi- 
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ation (14, 15). The (frequency-dependent) brightness temperature is given by: 



Here S v is the observed flux density at frequency v, 9 is the angular size of the 
source, and ks is Boltzmann's constant. At brightness temperatures above ^ 
10 12 K, multiple inverse Compton scattering of synchrotron radio photons to 7-ray 
and X-ray energies begin to rapidly dominate the luminosity. This is known as the 
"inverse Compton catastrophe", and effectively determines a maximal brightness 
temperature for a radio synchrotron source. 

Our most direct constraint comes from our VLBI detection on 2011 April 
1 (1), where we have measured S u {v = 8.5 GHz) = 1.7mJy. If we assume T B ^ 
10 12 K, we infer the angular size of the radio emitting region must be 9 ^ 3 /ias, 
or r = 9dx/2 ^ 3 x 10 16 cm. Assuming that the outflow began at the time of 
the initial Swift detection of Sw 1644+57 (201 1 March 25), this corresponds to a 
mean apparent velocity of (v) = r/At ^ 1.2c (ie., mild superluminal expansion). 
We note that this is entirely compatible with the fact that we have yet to (spatially) 
resolve the radio emission from Sw 1644+57: the upper limit on the angular size 
of 9 < lmas corresponds to an upper limit on the mean apparent velocity of 
(v) & 400c. 

For relativistic sources, Equation 1 must be modified to account for the dif- 
ference between properties measured in the rest frame of the emitting fluid and 
by an observer at rest with respect to the outflow. Accounting for the appropriate 
Doppler factors as described in ref. (15) (their Equation 3), our VLBI detection re- 
quires T = 1/ y/1 — (3 2 ^ 1.9; in other words, the outflow must be at least mildly 
relativistic. 

Separate from the above arguments based upon our VLBI measurement, we 
can use the observed variability of the radio counterpart to place constraints on 
the source size. The radio transient varies by approximately a factor of two over 
a time span of St ~ 7 hours in observations reported by the Owens Valley Radio 
Observatory [OVRO; (16)], and by ~ 50% over a time span of St fa 6 hours from 
data obtained at the Mullard Radio Astronomy Observatory [MRAO; (17)]. We 
conservatively adopt a variability timescale of 1 day for the radio (15 GHz) flux to 
change significantly. This variability could either be due to interstellar scintillation 
(ISS), as has been observed in the radio afterglows of gamma-ray bursts (18), or 
could be intrinsic to the source. 

In order to observe such strong variability due to scintillation, the observed ra- 
dio bands would need to be in the strong, diffractive regime, where the coherence 
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length scale is smaller than the Fresnel scale. This requires: 



--^3.4( IFn5 5M_) 6/,7 ^) 5/,7 GHz ,2) 

Here SM is the electron scattering measure, and D is the distance to the scat- 
tering screen. Using the electron distribution model of (19), we infer SM = 
1CT 3 - 66 m~ 20 / 3 kpc. For typical scattering screen distances (D fa 0.1-1 kpc), the 
scintillation frequency is thus z/ ss <^ 9 GHz. It seems likely, then, that the 15 GHz 
observations are not affected by strong ISS. However, we caution that the scatter- 
ing measure is only coarsely mapped, and small scale features can lead to large 
variations in SM. 1 

If the variability is indeed intrinsic to the source, we can use the observed 
timescale to independently constrain the source size. Based on causality argu- 
ments, the emitting region must be smaller than the light crossing time, and so 
9 ^ TcSt/dA ~ 0.2r/ias) at v = 15 GHz. Associating this angular size with a 
(relativistic) brightness temperature, we find T ^ 10, broadly consistent the limits 
on the angular sized and Lorentz factor inferred from the VLBI measurement 2 . 



1.4 Timescales and Accretion Rates for TDFs 

After the star is disrupted, a sizable fraction of its mass is placed onto highly 
eccentric orbits, which return it to the vicinity of the black hole on a timescale 
(22): 



^"fallback 




days, (3) 



where in the last expression we consider a solar-type star (R+ = R Q ) with a 
perinigricon distance 3 R p = 6GMbu/c 2 . The fall-back accretion rate peaks at 



'We note that there are no known supernova remnants (presumably arising in a particularly 
dense region of the ISM) in the vicinity of Sw 1644+57 (20). 

2 We caution that such causality arguments have difficulty explaining the observed fluctuations 
from the so-called "inter-day variables" (e.g., (27)) without resorting to extreme Lorentz factors 
(r ^> 100), and so may not necessarily be applicable to Sw 1644+57. 

3 We are considering here main sequence stars since giants will disrupt at larger radii and hence 
produce smaller accretion rates and (likely) less luminous jets. 
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t ~ tailback and declines oc t^ 5 / 3 thereafter (23, 24). Note that tf a iiback ~ days 
sets a characteristic timescale for the onset of bright emission, consistent with the 
observed onset of the activity associated with Sw 1644+57. 

The above estimates ignore the detailed physics of energy release near the 
MBH, assuming that the observed luminosity indeed tracks the instantaneous 
mass fallback rate. Even completely circularized accretion flows may be suscepti- 
ble to viscous and thermal radiation pressure-driven instabilities at high accretion 
rates (e.g., (25)). However, this is unlikely to be responsible for the short timescale 
variability, because the viscous timescale, t visc ~ a _1 (r /h) 2 Vt~ l is significantly 
longer than the dynamical time, unless both the disc viscosity parameter, a, 
is close to unity, and the disc geometrical aspect ratio is high, h/r ~ 1. A more 
likely source of short, dynamical timescale variability is the initial circularization 
of the disrupted material, which initially has high orbital eccentricity (26-28). The 
presence of the jet suggests that the BH may be rapidly spinning; thus there may 
be significant orbital precession near the horizon due to general relativistic effects. 
This precession could cause the falling-back stream of bound gas not to intersect 
itself on its first passage through perinigricon, instead taking several orbital pe- 
riods to shock and fully circularize its orbit (26); perhaps partial shocking and 
accretion at each pericenter passage gives rise to the observed X-ray variability. 

One plausible (though hardly unique) scenario that broadly describes the ob- 
served lightcurve is as follows: [1] t <; 10 5 s: initialization of a highly irregular 
jet; [2] t ~ 3 x 10 5 s: disc circularization phase, during which several large flares 
are powered as the bulk of the disrupted material on an eccentric orbit passes 
through the perinigricon (which probably strongly modulates the jet power and 
the supply of seed photons); [3] longer-term viscous evolution of the disk, char- 
acterized by a steadier consumption of material, fed by the continuing fallback of 
the disrupted stellar material, as described by equation 3; [4] a decrease in the ac- 
cretion rate from fallback material, leading to the dimming of the transient. Note 
that in this scenario the characteristic variability timescale is predicted to increase 
during the circularization phase. Circularization is probably mediated by strong 
shocks and may require up to several several local orbital times to complete. Thus, 
material closest to the SMBH circularizes first, while material further out requires 
a correspondingly longer time. This behavior is qualitatively consistent with the 
observed X-ray lightcurve. 
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1.5 Rate Estimates under the TDF Hypothesis 



We can obtain a model independent estimate of beaming by comparing the detec- 
tion rate of TDF events by Swift with the estimated occurrence rate. Assuming a 
tidal disruption rate per galaxy -Rtdf = -R-slCT 5 yr" 1 and adopting a local MBH 
density ~ 1(T 3 - 1CT 2 Mpc" 3 (e.g. (29)), the predicted TDF rate out to the red- 
shift z = 0.3534 (co-moving volume ~ 11 Gpc 3 ) is ~ 10 2 — 10 3 i?_ 5 yr" 1 . Since 
Swift has triggered on only one TDF in fa 6 years of continuous monitoring 4 , the 
required beaming fraction is fb ~ 2 x 1CT 3 — 2 x lO^^Rzl- Adopting typical val- 
ues for i?_ 5 ~ 1 — 10 predicted by theoretical studies (31, 32) and consistent with 
the TDF rate estimates from ROSAT (33), XMM-Newton (34), and GALEX (35) 
and constraints on off-axis radio events (36, 37), we find that f b ~ 10~ 4 — 10~ 2 . 
As discussed in the main text, this beaming is broadly consistent with the typical 
opening angles of blazar jets. This simple estimate above may need to be altered 
if the disruption rate or SMBH density is higher than our adopted fiducial range 
[as may be the case if the galaxy is small and compact and/or underwent a recent 
merger; (38, 39)]; or if only a fraction of on-axis TDF events produce nonthermal 
emission [as may be the case if another parameter other than the accretion — such 
as the SMBH spin — is necessary for jet production; e.g. (40)]. 

1.6 Argument for Separate Radio and X-Ray Emission Radii 

We now show that the radio emission and the X-rays (at least the variable contribu- 
tion) are unlikely to originate from the same spatial location. This lends credence 
to the suggestion in the main text that the X-rays are of "internal" origin, whereas 
the radio (and possibly IR) emission may originate from larger radii in the forward 
shock afterglow. 

If the emitting material has a Lorentz factor Tj, then from the minimum host- 
frame variability timescale t va r,min ~ 78 s, we infer the characteristic radius of the 
X-ray emitting region: 

Rx ~ cr 2 t var , min w 2.3 x 10 14 r 2 cm, (4) 

where Tj = 10r 10 . 

Separately, the radius of the radio-emitting region can be derived using the 
fact that any (incoherent) emission process cannot be brighter than the blackbody 

4 Note, however, that given the detection significance of the image trigger and the brighter later 
peaks in the light curve (30), this transient would have been eventually detected by Swift/BAT 
within a larger volume. 
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(BB) flux. Hence we demand that 

Fu ohs < F VobstBB = = ^^^ bs7cme (l + z f (5) 



167F „2„t,2„ 



L 

where 

Lu obs ,BB = ^r'^cT'-Uu, (6) 
is the specific BB luminosity and 

_ 87TV 2 kT e _ 87r^ 2 bs fcT e 

Uu ~ c 3 ~ c 3 rf () 

is the rest-frame BB energy density, v{y oh ^) are the rest-frame (observer) frequen- 
cies, and kT c fa 'j e m e c 2 is the "temperature" of the radiating electrons, corre- 
sponding to the characteristic (minimum) random Lorentz factor of the emitting 
electrons 7 e . 

If the radio emission is synchrotron, the peak frequency occurs at 

Vm — — , \o) 

2irm e c 

where B is the magnetic field strength in the emission region. If we assume that a 
fraction e B of the isotropic jet luminosity Lj is carried by Poynting flux, then the 
field strength at radius r is given by 



Combining equations (8) and (9) we solve for the electron Lorentz factor 



„ _ 10 2r-l/4 1/2 1/2 -1/4 (]Q) 
7e — J-U -£^47.5/15 i/ m,13 e B,-25 \ iKJ ) 

where Lj = 10 47 - 5 Lj i47 .5 ergs s" 1 , r = 10 15 r 15 cm, e B = lCT 2 e B ,-2> and u m = 
10 13 z/ m , 13 Hz. 

Now using equation (10) with (5) to constrain the emitting radius 

where i^ obs = i^ obsjmJy mJy, and z/ obs = lOz^o GHz. 

From the ~ 1.7 mJy VLBI detection at u w = 0.85, equation (11) shows that 
the emitting radius r^2x 10 16 cm, rather insensitive to the precise values of the 
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jet luminosity, synchrotron peak frequency, and ee- This constraint agrees with 
the minimum radius derived in §1.3, but remains valid even for large values of 
iy Note that the minimum radius of the radio -emitting plasma is incompatible 
with the X-ray radius Rx (eq. [4]) unless Tj ^ 100, significantly higher than 
values typically inferred for AGN jets. If the jet Lorentz factor were indeed this 
high, it would furthermore imply that the cooling timescale of the X-ray emitting 
electrons is much longer than t varim i n , a contradiction (41). We conclude that the 
radio emission and the X-rays (at least during the flaring state) originate from 
distinct radii. 

1.7 Origin of Inverse Compton Seed Photons in the Two Com- 
ponent Model 

One possible source of up-scattered photons in the two-component model is low- 
energy synchrotron emission itself (synchrotron self-Compton; SSC). In Figure 
1 of the main text we show model SEDs for the quiescent phase which include 
synchrotron-SSC emission. In general we find that pure synchrotron-SSC models, 
in which the source originates from sufficiently compact radii to explain the IR/X- 
ray variability, cannot by themselves explain the large ratio between the X-ray and 
IR luminosities. In addition, the peak flux of the second Compton peak is only 
marginally consistent with the Fermi upper limits. 

A more promising source for the X-ray emission (especially during the flaring 
phase) is external IC, as is thought to characterize the brightest blazars (e.g. (42)). 
Since this can in principle occur at smaller radii near the SMBH, it may also ac- 
commodate the flaring. What is the origin of the external photons? On timescales 
of just a few days, emission from the disk has not had sufficient time to establish 
a Broad Line Region (which is generally expected to form at r ~ 10 17 cm). Nev- 
ertheless other sources of soft photons exist within a distance r ~ 50r s ~ 10 14 
cm. One source is the thermal emission from the accretion disk, which is charac- 
terized by a temperature ~ 10 — 100 eV and luminosity ~ L c dd ~ xo 44-45 erg 
s -1 . However, for the jet to efficiently scatter the disk flux, then full acceleration 
to Tj ~ 10 must occur within a distance ~ 5 — 10r s , since otherwise the disk 
photons originate from the back and do not provide an efficient source (e.g., (43)). 

Another source of external photons is the shock where the bound gas falls back 
to perinigricon. While the fallback rate is super-Eddington, radiation pressure at 
the shock likely drives a powerful outflow of gas (11, 44) that can extend out to 
~ 50r s , a distance sufficient to extract sufficient energy from the jet via external 
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IC to power the observed X-rays. If the density at the shock is high enough to 
produce many photons, the radiation reaches thermal equilibrium and peaks in the 
optical/UV. IC scattering by even cold electrons in the jet will boost UV photons 
by a factor of at least r? ~ 100 (or higher if the scattering electrons are hot), 
resulting in hard X-ray emission. However, for M B h = 1O 7 M , the density in 
the shock may be too low for the radiation to reach thermal equilibrium; in this 
case, IC in the vicinity of the shock produce photons with energies ^ 10 — 100 s 
keV (44). These could then be further upscattered by the jet. We conclude that 
external IC is also a viable explanation for the X-ray emission. 

1.8 Forward Shock Emission Predictions 

Adopting the X-ray emission as a proxy for the jet power and duration, we as- 
sume that the bulk of the energy is ejected over a timescale T C j ~2 days, with 
average isotropic-equivalent luminosity Lj = 3 x 10 47 L 4 7. 5 erg s" 1 . The total 
energy release during this initial period is E- mo ~ 10 53 erg, after which the X-ray 
luminosity drops by a factor of several. 5 At a large distance, we model the jet as 
a homogeneous shell of thickness 5 = cT C j = 5i 6 10 16 cm, which propagates into 
an external medium of density n cxt (r). The circumnuclear gas is expected to be 
rather dense out to the Bondi radius R h due to e.g. gas accretion (e.g. (45)). For 
simplicity, we assume that the density scales oc r~ 2 for r < R h = 0.1 pc, but then 
smoothly connects onto a constant density profile n cxt = I0n w cm~ 3 further out. 
Our results depend rather weakly on the details of the profile as long as n ext ~ 10 
cm -3 for r ^ 0.1 pc. 

The interaction between the jet and the external medium results in a charac- 
teristic forward/reverse shock structure. The Lorentz factor of the shocked fluid 
is (46) 

r sh =(^) 1/4 = 3.< 4 5 <o 1/4 , (12) 

where rij = Lj /4irr 2 T? m p c 3 is the number density in the rest frame of the jet. 
Note that in the last expression we have assumed that r < i?b, in which case T s h 
is independent of distance because both the density of the jet and external medium 
scale oc r~ 2 . During this stage the jet is coasting, such that the emission from the 
shock is observed at a time t Q bs ~ r /^ 2 h c °z r ' beginning after the onset of the jet 
activity. 

5 Note that jet activity persists by injecting more energy at a lower rate, which may contribute 
to the late-time light curve. 
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In modeling the emission from the external shocks, shown in Figure S2, we 
adopt the gamma-ray burst afterglow model of Sari et al. (47). We assume that a 
fraction e c and e B of the jet energy is placed into that of non-thermal electrons (that 
follow a power-law distribution with p — 2.5) and the magnetic field, respectively. 
The time dependence of the self-absorption, characteristic, and cooling frequen- 
cies of the forward shock emission are, respectively, given by u s& oc u m oc f -1 , 
and v c oc t. The flux at the characteristic frequency v m is constant with time at 

1/2 1/2 

the value F vm = 100e B _ 2 £ , 53 n 1 mJy. For this choice of parameters, emission 
from the reverse shock is subdominant to the reverse shock because (a) the rel- 
ative Lorentz factor between the the shocked fluid and unshocked jet is modest 
r re i = 0.5(Tj/T s h + r sh /Tj) ~ 1.6; and (b) the reverse shock emission suffers 
substantial self-absorption. 6 

The reverse shock crosses the ejecta at a distance R CT ~ r 2 h A ~ 2 x 10 17 L l J 7 2 b n 
cm (46) similar to the assumed Bondi radius R h . At radii r > R cr , the blastwave 
propagates into a constant density medium and reaches a self similar evolution 
T s h oc r -3 / 2 . This transition occurs at an observer times that is slightly longer than 
the duration of the jet activity [defined as when most of the energy is ejected; (49)]. 

Once the self-similar evolution is established, the evolution of the synchrotron 
spectrum is described as in Sari et al. (47). When the synchrotron frequencies 
^ sa < v m < v c they find that u sa = 2 x 10 10 e B ( 5 .2^-i-^53 n i 5 Hz > u m = 
3 x l^e^el^EliX^, and , c = 8x lO^e-^^^nr 1 ^ 2 . These ex- 
pressions are valid for observer time t ^ cTj ~ 2 days. The flux at the charac- 
teristic frequency is F Vm = nOe^^E^n 1 / 2 ^^ mJy, where we have applied a 
correction (^r s h) 2 ^ 1 to the Sari et al. (47) expression to account for the fact the 
observer already observes the edges of the jet because (unlike GRB afterglows) 
we have assumed that initially the Lorentz factor and opening angle are related by 

e 3 = i/ry 

6 Note that Giannios & Metzger {48) assumed a much longer event duration, i.e. months not 
days. As a result they concluded that the reverse shock was relativistic and dominated the emission. 
This limit may still apply to Sw 1644+57 at late times if a lower luminosity jet persists. 
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Fig. Sla (top): Variability analysis of the X-ray light curve (gray with green 
spline overlay) observed with the Swift XRT. We define the variability timescale 
t var as the time between count rate rising (red) or declining (blue) by a factor 
of 2. Upper limits on i var can occur between gaps in the observing coverage 
by the satellite, (bottom) Cumulative rate of significant rising and declining 
events from t + 10 3 sec to t + 1.4 x 10 6 sec. A number of the measured 
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100 - few x 10 s ) occur in the timespan from t + 3 x 10 4 sec to 



t + 3 x 10 5 sec. 
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Fig. Sib: Inferred power density spectrum of the Swift/XRT light curve. The 
solid straight line shows the fit to a single power ("red noise") along with the 
1 a uncertainty range (dotted lines). 
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Fig. S2: The predicted light curves from the forward shock synchrotron 
model at typical radio v — 1 GHz, millimeter v — 100 GHz, and IR fre- 
quencies v = 10 14 Hz for E 53 = 0.3. The observed radio and IR light curves 
are overplotted. A prediction of this model is that the 100 GHz/IR emission 
will decline at late times. Note, however, that if the X-ray emission from 
J164449.3+573451 continues at a similar level, soon the energy injected by 
the quiescent phase will eclipse that injected in the first few days. This addi- 
tional energy is not accounted for in this model and may cause the light curves 
shown to decline less rapidly at late times. On timescales ~months-years the 
light curve will revert to the predictions of Giannios & Metzger (48). 
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Fig. S3: Two purely-synchrotron models fit to the instantaneous SED at t + 
2.9 day, following from Fig 2. of the main text. The dark blue line shows a 
simple synchrotron spectrum matched to the data. A large extinction column 
(A v > 10 mag) is needed to explain the relatively faint optical/NIR fluxes, 
and the fit is poor. In addition, the large value of the cooling frequency v c in 
this model requires a source size inconsistent with the observed variability. 
The red line shows a two-component synchrotron model, consisting of (1) a 
more compact emitting region more consistent with the optical/NIR photom- 
etry and variability timescale and (2) a more radially-extended synchrotron 
component, which represents synchrotron emission from the afterglow (see 
Section 1.8 and Fig. 1 of the main text). 
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